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: Collision Trajectories in the 


Three-body Problem 


E. T. Benedikt* 


Abstract 


The connection between the analytical regularity of the 
trajectory of a missile moving under the action of two at- 
tracting masses (in practice the Earth and Moon, or the Sun 
and a planet) and the collision between the missile and the 
center of one of the bodies is noted. By applying a trans- 
formation due to T. Levi-Civita, it is possible to regularize 
collision trajectories, and to obtain conditions which the po- 
sition and velocity of the missile must satisfy (at a point 
sufficiently near an attracting body) in order that the tra- 
jectory of the missile pass through the center of the body. 


Introduction 


The determination of collision trajectories, that. is, 
of trajectories passing through the centers of celestial 
target bodies, can be obviously regarded as one of the 
basic problems of Astroballistics. t 

Results of astronautical interest obtained so far in 
this field appear to be of either elementary qualitative 
nature, or brute-force automatic numerical computa- 
tion. Analytical results would be much more intellectu- 
ally satisfactory as well as useful in connection with 
such problems as the choice of trajectories, assessment 
of accuracy of ballistic missiles, etc. It seemed there- 
fore in order to draw attention to certain very funda- 
mental results obtained over 50 years ago by the great 
mathematician, T. Levi-Civita. Possibly due to lack of 
a suitable field of application, his results were not 
further used (to the knowledge of the writer) except as 
dynamical theorems pertaining to the three-body 
problem. It appears, however, that further develop- 
ment and possible extension of Levi-Civita’s methods 
might well yield analytical results of great usefulness to 
astronautics. 


Analytical Singularity and Collision Trajectories 


The three-body problem, i.e., the integration of the 
equations describing the motion of three particles under 
the action of their mutual gravitational attraction con- 
stitutes, as is well-known, one of the classical problems 
of celestial mechanics. It is equally well-known that the 
above system of differential equations cannot be inte- 
grated by quadratures. As long as the trajectories which 
are being sought do not pass through one of the centers 
of attraction, i.e., outside the domains in which the 
differential equations of motion are affected by singulari- 

* Group Physicist, Missiles Research Section, Douglas 
Aircraft Company, Inc. 

_ + We choose to designate with this name the engineering 
of unguided ballistic interplanetary missiles. 


ties, it is always possible to obtain a sufficiently approxi- 
mate particular integral in the form of a power series of 
the variables, by the usual procedures of Cauchy. These 
procedures fail for trajectories passing through one of 
the attracting centers, the failure being actually limited 
to a neighborhood of this body (the size of this neighbor- 
hood depending, of course, upon the desired approxi- 
mation). 

A connection (first noted and investigated by Pain- 
levé', pages 582-86, exists therefore between the possi- 
bility of (binary or general) collisions in the past or 
future history of a system of three gravitating particles, 
in a given initial configuration and with given initial 
velocities, and the possibility of obtaining integrals in 
the form of holomorphic functions of the ensuing 
motion. The possibility of obtaining integrals of col- 
lision trajectories (regarded as singular trajectories) 
was examined more closely in 1903 by Levi-civita’ for 
the case of the plane reduced three-body problem. (The 
term “reduced” implies that one of the particles of the 
system is a planetoid, i.e., of mass so small as not to 
influence the motion of the remaining bodies. An 
analytical condition permitting identification of the 
collision trajectories emerged from this analysis. Closely 
connected with this work is the subsequent notable 
analytical results obtained by the same author on the 
regularization of the two-dimensional restricted three- 
body problem’. This “regularization” consists essen- 
tially of a change of independent and dependent vari- 
ables such that the resulting equations remain regular 
even at points occupied by one of the attracting centers. 
A very important feature of Levi-Civita’s method of 
regularization is that the transformed equations pre- 
serve the canonical form typical of the equations of 
dynamics. Subsequently (1912) Sundman’ succeeded 
in regularizing the three-body problem in its most 
general (three-dimensional) form. For this result, 
Sundman was awarded in 1913 the Pontécoulant prize 
of the Paris Academy of Science, and hailed by some 
mathematicians (and even by the press) as the dis- 
coverer of the solution of the three-body problem. 
However, his methods are very lengthy and inelegant, 
and the regularized equations established by him suffer 
from the drawback of not being representable (as Levi- 
Civita’s equation ) in the canonical form, the theoretical 
results and methods of analytical dynamics being there- 
fore inapplicable for further developments. 

These developments, as well as the later contribu- 
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tions of T. Levi-Civita” ° in the field, are however, not 
of direct importance in astronautics, where interest is 
obviously limited to the restricted problem of ballistic 
space ships functioning, of course, as planetoids. 

We shall proceed to rederive the basic results ob- 
tained by Levi-Civita in a manner following the formal 
methods of dynamics. { 


Lagrangian Function for Relative Motion of 
Planetoid 


Let us consider two masses M, M, assumed concen- 
trated in two points 0, 0; and subject exclusively to 
their mutual Newtonian attraction. If further the 
distance 00; = a is constant, the absolute motion of the 
two masses must, as is well-known, be the resultant of a 
uniform rectilinear motion of their barycenter 0 placed 
on the line 00; at a distance 


7 = M,a/(M + M;) (1) 


from 0; and a uniform rotation about the barycenter 
with angular velocity 


my = Vk(M + M,)/a (2) 


about an axis perpendicular to the line 00; , k being the 
gravitational constant. 

Consider now a third body (planetoid) moving under 
the action of the gravitational field of the two attracting 
centers. We shall restrict our considerations to the plane 
in which these centers move. The Lagrangian function 
of the planetoid will then be given by the difference 


Lo ACU 


(3) 
between its kinetic energy per unit mass 3 = (4)u9° 
(where vp is the absolute velocity of the planetoid) and 
its potential gravitational energy 


U = —kM/r — kM,/n (4) 


similarly per unit mass, 7, 7, denoting the distances be- 
tween the planetoid and the attracting centers 0, 0, . 

Introducing a system of polar coordinates 74 having 
origin at the barycenter 0 of the system and in which 
the angles are reckoned from the line 00, (lig. 1) we 
shall have 

Si (4) [ro aie ro (00 fe no)’| 

Expressing this function in terms of a new system of 
polar coordinates r, 6 having center in 0, and in which 
the angles are reckoned from the line 00,, the above 
expression can be written in the form 


2 22 P 2, 2, 2 pes 
3 = (s)[F + rf + 2m 16 + m (rr — 2rF cos 6)] 
.2 ar 2 Oke 
= ()[F + 7(@+ mm) — 2norF cos 6] 
where inessential constant and linear (with respect to 
the velocities) terms have been disregarded. Substitut- 


{t A treatment following somewhat more closely the original 
analytical developments of Levi-Civita, contained in the 
(not easily accessible) reference 2, can be found in reference 7. 
A discussion of Levi-Civita’s regularization procedure, par- 
ticularly from its geometrical aspects, is given in reference 8. 
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Fig. 1. Absolute and Relative Motion of Planetoid 


ing this result together with (4) into (3) dividing 
through by ma and employing identity (2), as well 
as the nondimensional coordinates and parameters de- 
fined by the substitutions 


Fw IGE 
p= Ta 
A = r,/a 


the Lagrangian of the planetoid can be expressed in the 
form 


£ = (4)[o" + p (6 + 1)] —§ (5) 


where 
F(p, 0) = upcos@ + U 

= up cos 8 — v/p — w/A 
Here and in what follows, primes will denote derivation 
with respect to the parameter 7. In (6) we have further 
put »y = M/(M + M,) =1— 4. 
Hamiltonian Function of the Planetoid 

From (5) we obtain the expressions 
Pp = OL/dp’ = p’ 
Po = 02/36" = p'(6' + 1) 


(7) 
(7") 
for the momenta conjugate to the polar coordinates. 
From (7’) it follows that the angular velocity 6’ can be 
expressed in terms of the momentum py by means of the 
relation 
: 

O > piy Dares (8) 
By means of (7) and (8) it is possible to express the 
Lagrangian (5) in the form 


£ = (3)[pe + (1/0')p'] — U (9) 
The Hamiltonian function defined by 
= ppp + po — & (10) 


with the velocity expressed in terms of the momenta 


will on account of (7), (8), and (9) be given by 


K = [{p, + pl[(1/p)p — 1}} + UI] 


(11) 


j 
: 


(6) 


where 
UO =F — ($)p = Volr) + wVi(r, 6) (11’) 
The Jacobian Integral 


Since the Lagrangian (8) does not depend explicitly 
upon the independent variable, the expression 


p'dL/dp' + 0'dL/a0' — & 

must remain constant and thus represent a (first) 
_ integral of the equations of motion. Indicating with 
w= Vp? + p76 (12) 


the (non-dimensionalized) speed of the planetoid (con- 
nected to its actual speed v by the relation v = m aw) 
we obtain the so-called Jacobian integral 


Qw+o= K (13) 


of the reduced three-body problem. The existence of 
this integral follows even more directly from the Hamil- 
tonian point of view, where the fact that the inde- 
pendent variable does not appear in it ensures the 
constancy of this function. The Jacobian integral can 
thus be expressed in terms of the momenta in the form 


KH = (4){pe + rl(1/p)pe — 1p} +V=K (14) 


_ The left hand members of relations (13) and (14) and 
consequently the constants appearing in the second 
member are numerically identical. From (13) and 
(11’) we obtain immediately the formula 


w= V/2{K + (3)p? + »/p + w[(1/A) — pcos6)]} (15) 


Transformation of the Radial Coordinate 


Replacing the radius vector p by the new variable o 
defined by 


g=|Ve| p=o (16) 
the Lagrangian function (5) becomes 
& = (4)[4o'o” + o°(" + 1)'] —F(c, 6) (17) 
The momenta conjugate to o, @ will become 
De = OL/d0' = 40'0' (18) 
po = 0/00’ = o (6 + 1) (18’) 


the velocity components being conversely expressible in 
terms of these momenta by means of the relations 


o’ = pa/4o° (19) 
6 = 9 0/0 — 1 (19’) 
The Hamiltonian function will consequently assume 
the form 
5 = poo + ped’ — £& 
{(40") pe’ + [(1/o*)pe — oF} + 0 
where U is defined as before by (11’) which now 
becomes 
U = —v/a’ = (4)o* — pl(1/A) — o° cos 6] (20’) 


(20) 


The Jacobi integral will be in this case 

K = (3){(40°)pe + [(1/o") pe — 0°} (21) 
— r/o — p/A + po’ cos 6 — (4)o* 

Regularization of the Equations of Motion 


The motion of the planetoid will be defined by the 


system of differential equations 
o' = 05/dp, 6’ = 93C/dpe 
(22) 
De = —O5/dc py’ = —O5C/98 


where 5: = 3C(c, 8; p., po) is the Hamiltonian (20). 
Introducing the new dependent variable defined by 


ds = dr/o (23) 
equations (22) assume the aspect 
do/ds = o 03C/dpe do/ds = o A3/dpe (24) 
dp,/ds = —o dK/da dps/ds = —o d35C/d0 


Due to the existence of the Jacobian integral (14) 
we have 


d(a 3) /de = 5da'/do 
d(a 5C)/do — 25C 
d[o' (Ke — K)]/do 


Since o is independent of co, ps and ps we can write 
furthermore 


o O5C/d0 


o 05/00 = d[o (Ke — K)]/d0 
¢ 05/dp. = Alo (KH — K)\/dpe 
o 05C/dp = do (KC — K)]/dpe 


‘The function 
2 
=¢(K — K) 
: ; 1 2 3)2 (25) 
= (5) (2) pe ([Gfe) pr ro ire 
where 
@ = +v+ Ko’ + (2)0° 4+ ul(o’/A) — o cos 6] (25’) 
can be thus taken as the new Hamiltonian function 


by means of which the motion of the system is defined 
by the canonical system of equations 


da/ds = 0G/0p. d6/ds 
dp,/ds = —0G/0. dpe/ds 


It should be noted that the function G defined by (25) 
is devoid of singularities at the origin (ep = « = 0) 
of the coordinate system. 


I 


0G/dpe (26) 
—ag/ae 


Hamilton-Jacobi Equation 
Putting 


OW /do 
aw /d0 


me (27) 


Pe 


into the regularized Hamiltonian and equating the 
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result to a constant @ (which, in our case, is numerically 
zero) we obtain the Hamilton-Jacobi equation 


(4){4(9W/de)’ + (8W/0d8 — o°)’} 
=G+v+Ko + (4)o°+u(o0/A — oc cos8) (28) 
A complete integral of this equation will be of the form 
W = Wa, 0, G, h) 


h being a constant of integration. According to the 
well-known results of the Hamilton-Jacobi theory, the 
orbit will be defined by the (implicit) equation 


(OW/dh) a0 = % (29) 


where @ is an additional constant (determined by the 
initial conditions), whereas the kinematics of the 
motion will be given by 


Ss — SS = (OW /dG) go (30) 


It is apparent from the above results that if we set 
G = 0 in (28), formula (29) will remain unaffected; 
equation (30) must however be replaced by the equiva- 
lent expression 


S SS Gt) = (0/dv)(W) g=0 (30’) 


Determination of the Characteristic Function 


From equation (28) we obtain 
IW /do0 = F+V/8[P — (4)(0W/cd0 — o*)2| (31) 


A complete integral of equation (28) can be obtained 
by the classical methods of Lagrange, Monge, and 
Charpit* by considering an additional relation of 
the form 


dw /d0 = O(a, 0) ED) 


the function © being subjected to the condition of 
integrability 


(d/d0) (dW /d0) — (d/00)(0W/dc) 
= 00/00 F(2/cR){[(O/a) — «°|(00/00) (32) 
— o0°/00} = 0 
where we have put 
R= | V2 — ($)[(O/c) — oF]f*| (32’) 


The characteristic function will then be given by the 
integral 


Aes il (+28 do + © dé) (33) 
extended over an arbitrary path of integration. 


Kinematical Result 


If an integral of the form 


o(c, 0,0) =h (34) 
* Equation (28) cannot be integrated by the usual methods 
of separation of variables except in the case y = 0 cor- 


responding to the two-body problem. 
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; : f 
of the quasi-linear partial differential equation (32)8 


is available, the equation of the orbit of the planetoid | 
will be given, according to equation (29), by i 
aw/ah = / {+ (2/8)[(@/a) — 0°(00/ah) do (a5) 
If the integral (34) appears in the more special form 

(oc, 6,0) = O — O(a, #) sh =0 (34% | 


we have 00/dh = (06/80) = —1 and consequently _ 
equation (35) becomes 


at he + | — (2/oR)[(@/s) — o*)| do (36) 8 


The equation of the orbit can also be expressed in the — 
differential form | 


d0/do = ¥(2/oR)[(@/c) — ol (37). 


The angle y formed by the velocity of the planetoid and 
the radius vector satisfies the relation tan y = pd6/dp 
which, applying transformation (16) becomes 


tan y = (4)ed0/do 
On account of (37) we have 
tan y = #(1/@)(O/o° — o’) 


Still assuming an integral of (32) expressible in the— 
form (34’) we obtain from (30’) and (33) the equation 


8 — 8 = #(2/R){1 — [(O/c) — o']8O/cdv}do 
+ (d0/dv)dé@ 


This expression can also be put in the form 


peu [ + ((2/R){1 — [(@/o) — o°](a@/oar)} 
+ (d0/dv)(dé/do) do 


where the path of integration must now be restricted 
to an actual are @ of trajectory. Using (37) we have 
therefore 


S—% = = 2/R do (38) 
: 
or, in differential form 
ds — #F 2da/R 


Since, on account of (23) 


2 
ron= | ods 


it follows that the time dependence of the motion is 
defined by the equation 


r= 1 = ¥[ 20°/8 do (39) 


Differential Equations of Trajectory 

The preceding developments depend upon the 
knowledge of an integral of equation (32). If this equa- 
tion were to be solved by the conventional method of 


characteristics, it appears immediately that one of the 
(ordinary) differential equations for these lines, namely 
the one defining d6@/do coincides with equation (37) 


whereas the remaining equation, that is 


dO/do = +2(0°/00)/+/2 — [(@/c) — a]? (40) 


can be regarded as auxiliary equation defining 0*. 


determined directly from equation (32) 
manner, of course, as to fulfill the same initial condi- 


_ Equations (37) and (40) constitute the system of 


differential equations, an integral of which—particu- 
larized for given initial conditions—would coincide 
with the trajectory defined by equation (36), the 
function © appearing in this equation having been 
(in such 


tions). Therefore, application of the method of char- 


acteristics to the present problem would be essentially 
equivalent to the direct integration of the system of 
“canonical equations (27).7 It is therefore apparent 


that the methods of integration described under the 
section on “Determination of the Characteristic Func- 
tion” is valuable only as long as a direct determination 
of the function © is possible. This is fortunately the 
case, at least for trajectories through the attracting 
center. 


Collision Conditions 


From (40) it follows that if the quantity of dO0/do 
is to be real for ¢ = 0, the limiting value of (Q/c) must 


be finite, from which it follows that the condition 
C0) ean) (41) 


is necessary if the trajectory of the planetoid is to pass 
through the attracting center. Since 


A =~V/1 — 2pcos6 + p? = V/1 — 202 cos6 + o4 (42) 


we obtain from (25) 


00/86 = —po[1/A® — 1] sin@ = oQ (48) 
where 
3 © is ie 2 7 
= —(3)p[2 cos 6 — (1 — 5 cos 8)o 
Q (ul (43’) 
+ 0 (4)] sin 0 
Inasmuch as @(0) = v and Q(0) = —3p cos @ are 


finite, it follows from (40), (41), and (43) that we 
must have 


(dO/do)o-0 = 0(6) (44) 


This suggests the possibility of fulfilling equation (32) 
by assuming 


© = a'W(a, 0) (45) 


* The characteristics of equation (32) can be represented 
in the space obtained adding to the plane of the motion a 


third dimension representing 6. The projection of the above 


characteristic lines, defined by equations (37) and (40) upon 
_ the plane of the motion define then the orbits of the planetoid. 


+ Equations of the trajectory obtained by elimination of 


the time from the canonical equations have been indeed 


systematically used by Levi-Civita! to establish the collision 


conditions. 


y being a holomorphic function of o, accordingly repre- 
sentable in the form 


v= Dal A)o" (45") 


Substituting (45) into (32) making use of (43) we 
obtain 
ody/do + TW = +2/R{o'lo*'y — 1lay/a9 — ©} (46) 


the factor o° common to all members having been 
eliminated. 
Since 


1/A — o cos = 14+ 0(4) 


we obtain from (32’) and (45) the expansions 


ns 2v + 2Ko° + wo [1/A — o cos 8] 
eh 1 fy/ — ol(o4y — 1)? — 1] 


1/+/2r{1 — [(K + p)/2r]o” + 0(4)} 


From this result and (43’) we obtain further 


(47) 


I 


Q/R = (8u/2v/2r) sin 6(2 cos 0 
—{1+ [(K + p»)/2y] cos 6 — 5 cos’ 6}o° + 0(4)) 


with the aid of these results and using (45’) we obtain 
the expansion 


+ (2/R){o'lo'y — llay/a9 — Q} 
= +(1/+/27) [By sin 6(2 cos 0 
—{1 + [(K + 2)/»] cos 6 — 5 cos’ Ofc’) 
 20°dWo/d0 + 0(4)] 
Using (45’) we have for the left hand member of (46) 


(48) 


oO 


cdw/do + 7W = >> (n+ 7)yno” 


n=0 


(49) 
Comparison of the coefficients of similar power of o 
appearing in (48) and (49) yields the determinations 

Yo = (4) (u/v/2y) sin 26 

yi = 0 

vo = + (u/12/2v){ — sin 0 

— 2[(K + yw)/r| sin 20 + 5 sin 36} 
v3 = (3u/35v) cos 20 


(50) 


for the coefficients of an expansion of y valid up to 
the third power of oc. From these results, (45’) and 
(45), we obtain the condition 


O = +(u/vV/2r)0°((#) sin 20 + (4 
{sin 6 — 2[(K + y)/r]sin 20 + 5 sin 30}0° 
= (3u/35v) cos 260°) + 0(10) 


to which the angular momentum of the planetoid must 
be subjected in order that the trajectory of the latter 
pass through the attracting center. 


(51) 
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Fig. 2-a. Required Transversal Velocity for Collision Tra- 
jectory of Planetoid (K = —.47895). 


Collision and Ejection Trajectories—Symmetry 
Relations 


Let us assume that the function © contained in 
(38) fulfills condition (51). The quantity 


do/dt = -Q@/2c 


will increase or decrease, that is, the planetoid will 
move away from or towards the attracting center 
according as to whether the sign affecting the radical 
® appearing in the preceding development is negative 
or positive. Therefore, at least in a neighborhood of 
the origin in which there are no reversals in the direc- 
tion of the motion, the above formulae will pertain to 
collision or ejection trajectories according as to whether 
@® is taken positively or negatively. In particular, ex- 
pansion (51) will define the condition for collision if 
taken with the upper set of signs, and ejection if taken 
with the lower set of signs. 

From (5) and (6), it is apparent that the Lagrangian 
function of the planetoid remains unchanged, if the 
variables —6@ and —7, are substituted for 6 and +r. 
This symmetry implies that a reflection about the line 
00, will change a collision into an ejection trajectory 
and vice versa. The same conclusions can be reached 
by considering that a change of @ into —@ in formula 
(51) will interchange the upper with the lower row of 
signs and thus collision into an ejection condition and 
conversely. On the other hand, formula (37) can be 
made to define a collision or ejection trajectory by 
substituting in it the appropriate expression for © and 
®, taking ® with the appropriate sign. Transition from 
one type of motion to the other implies thus a change 
in sign of the right hand member of equation (37) 
as well as of dé/dc, and therefore, of © in accordance 
with the above results. 

Finally, it is worth mentioning explicitly that the 
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Fra. 2-b. Required Transversal Velocity for Collision Tra- 
jectory of Planetoid (K = +.46058). 


role of 0 and 0; as respectively attracting and perturb- 
ing center are interchangeable. Results obtained for 
trajectories through 0 can thus be applied in trajectories 
through 0; by an obvious change of the meaning of the. 
various symbols. 


Alternate Form of the Collision Conditions 


From (31’), (45), and (19’), it appears that the 
angular velocity of the planetoid (relative to the 
synodic system here employed) will be given by 

= —-14+ O/o* = -1+4+ op (52) 
Making use of (50) and reverting by means of (16) to: 
the original coordinate p, it follows that in order to: 
collide (or to have collided in the past) with the at- 
tracting center 0, the rate transversal component. 
Wes = 70 of the relative velocity of the planetoid must. 
satisfy the condition 


we = —p + w/r/2vp"*((F) sin 26 
+7's{sin 6 — 2[(K + p)/v] sin 26 + 5 sin 36}p (53) 
F (34/35 v) cos 26p"”) 


where the upper signs refer, as before, to collision and. 
the lower to ejection trajectories. 


Analytic Representation of Collision Trajec- 
tories 

A graphical representation of the right-hand side 
of equation (53) is shown in Fig. 2-a and 2-b for the 
cases (of practical interest circumlunar ballistics) for 
which K assumes, respectively, the values —.47895. 
and +.46058, : 

On account of (32’), (25), (45), and (45’) we have 
the expansion 


1/R = (1/v/2r) a Cro” + 0(9) 


Cy =A 
Mec (Kn) 2p 
C2 = +(3)[(K + 2)/of? 


= 


= —{(u/8v)(1 + 3 cos 26) 
+(48)[(K + u)/v]"} 
Cs = 3u/16r{(K + w)/v — cos 6 
+[3(K + p)/r] cos 26 — (8) cos 30} 


ubstituting the above results into (36) and (39) and 
mtegrating, we obtain the approximations 


++/2/v be o” 09) | (54) 


6 — 6 


: sila sae V2» D Bo" + a9) | (55) 


efining trajectories through 0 in the neighborhood of 
shis point with 


ay = 3 

a, = 0 

as = —(K + p)/10» 

as = — (te) (u/-V/ 2) sin 26 (55!) 


ar = +(e5)[(K + u)/o 
a = F(4)(u//2r){— (7s) sin 0 
—+4((K + pw)/vr| sin 26 + +5 sin 36} 
nd 


Bs = a3, Bs = a5, 87 = a7, 841 = Be = Bg = 0 (55”) 


ollision Condition for Keplerian Motion 


For the case in which the perturbing mass is zero, 
he above results describe a Keplerian motion with 
respect to our synodic system of coordinates. In this 
ase, the collision condition (51) reduces to 


@.= 0 (96) 


ecording to (52’), this implies that in this case the 
ngular velocity of the planetoid about 0 must be —1 
elative to the synodic system of reference. Inasmuch 
the latter rotates with respect to an inertial system 
f reference with an angular velocity which, with the 
nits adopted here, is +1, it follows that in order to 
ollide, or to have collided, with 0, the absolute angular 
elocity of a planetoid must be zero. This means that, 
as known a priori, the unperturbed collision trajectory 
relative to an inertial reference system must be straight 
lines through 0. As apparent from (52) the absolute 
angular momentum 0 = o (6 + 1) = p (0 + 1) of 
such trajectories is zero. If » = O(v = 1), (25’) re- 
duces to 

@ 22 Ko a (h)e" (57) 


If moreover the collision condition (56) is fulfilled, 
(32’) becomes 


R= | V2 + Ko? | (58) 
From (56), (37), and (38) it follows also that 
dé/dr = 1, that is 
et ee 


The motion will then, according to (36) be defined in 
this case by 


ron = £2 | (odo)/VI+ Ke (59) 
or, on account of (16) 
rm = F(1/V2) | (Wo do/VT + Ke) (50')* 
Taking the integration constant ro as 0 (i.e., for a 
planetoid located at 0 for 7 = 0), we obtain thus 
(F1/v/2[(o/K) V1 + Ko? + (1/—K)*? 
sin +1/— Ko] for Ka) 
ee Oe for K = 0 (60) 
F1(//2)[(c/K) /1 # Ko? — (1/K)” 
sinh | 1/Ko] for kK > 0 


Trajectories of this type relative to a synodic system 
of reference are shown in Fig. 3 for typical cases of 
elliptic (K < 0), parabolic (K = 0), and hyperbolic 
(K > 0) trajectories.| The expansion in series of o 
of (60) coincides, of course, with the one defined by 
(55) and (50') for n — 0. 


Application to Lunar Ballistics 


The practical applicability of the foregoing develop- 
ments was tested on the case of cireumlunar trajectories. 

The calculations were made on a model of the Earth- 
Moon system in which the motion of the Moon around 
the Earth was assumed to be circular, all the perturba- 
tion due to other celestial bodies regarded as non- 
existent, and the value of uw was taken as 0.98787. 
Trajectories with initial conditions fulfilling equation 
(53) and for which the Jacobi integral assumed the 
values K = —.47895 and K = +.46058 were com- 
puted by means of the Runge-Kutta method on a 
type 704 IBM digital analyzer. Initial distances from 
the Moon of the order of six terrestrial radii were used. 

* Equation (59’) can be obtained immediately by separa- 
tion of variables from the principle of conservation of energy, 
which expressed with respect to an inertial system in the case 
at hand is given by the relation 

(4)p — 1/p = (4)p0? — 1/p0 

popo defining some initial position and velocity. It should be 
noted that the right hand member coincides with the con- 
stant K defined by (13) only if the angular velocity 0’ (rela- 
tive to the synodic system) is +1. 

+ The trajectory of the planetoid does not reach, is tangent 
to, or crosses the orbit of the (massless) satellite according as 
to whether A < —1, K = —0,orK > -1. 
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Fie. 4-a. Earth-Moon Collision Trajectory With Respect 
to Synodic System of Reference (K = —.47895). 


Points of departure located at greater distances seemed 
to fall—at least for the particular values of K employed 
here—outside the domain of convergence of the series 
(51) or (53).* Trajectories obtained by the procedure 
described above, shown as Fig. 4-a, 4-b, 5-a, and 5-b, 
pass extremely near the center of the Moon. In evaluat- 
ing the results, it should be kept in mind that the 
initial conditions were computed on the basis of an 
approximate analytic representation of the collision 
conditions. (The errors inherent in the numerical pro- 
cedure, estimated to be a fraction of a mile, were of a 
much smaller order of magnitude.) These findings 
seem to be quite encouraging, since they indicate that 
in spite of the large value of the perturbing mass (of 
the Earth), collision conditions can be adequately 
approximated within a domain of significant linear 
dimensions. 


References 
(See page 30) 
* For additional information regarding the question of con- 


vergence of the expansion of the collision conditions, see 
reference 7, 
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A Review of the Current Status of 
Satellite Attitude Control’ 


Robert E. Roberson} 


‘Abstract 


The literature pertaining to each of the major aspects of 
satellite attitude control is summarized and discussed. These 
include general descriptions, first principles and analytical 
formulations, control methods, sensors, actuators, synthesis 


} and mechanization, and performance analysis. The status of 


each area is examined and some problems are indicated where 
further work is especially desirable. 


Introduction 


Attitude control of a satellite vehicle is a problem that 
‘has not received much attention until relatively re- 
cently. Nevertheless, with the growing realization of its 
importance to the success of most satellite missions, a 
substantial body of literature on attitude control is now 
~ available. 

The purpose of the present paper is to summarize and 
‘discuss this literature. The framework for discussion is 
the set of major problem areas which characterize 
satellite attitude control. Contributions to each prob- 
lem area are listed, and the status of each problem is 
examined. Some areas are indicated where further work 
is especially desirable. 

The literature cited is of three basic types. The bulk 
of it specifically concerns the attitude control problem 
per se, or at least includes a significant amount of atten- 
tion to this as one area among several. Usually there is 
some brief mention of attitude control requirements in 
relation to the mission or application, but this is in 
passing and is in the nature of background rather than 
detail. On the other hand, other papers primarily con- 
cern the application, with passing mention of the need 
for attitude control. Only a few references of this kind 
are given, specifically (25), (26) and (27); undoubtedly 
a number of others could be discovered. One of these, 
(26), amplifies the discussion of attitude control, giving 
a brief description of a suitable attitude control tech- 
nique. The third category is represented by just one 
entry, (4). It is the complete discussion of the attitude 
control problem—requirements’ and system—for a 
specific mission; in this case, for a lunar instrument 
carrier. Similar studies do exist for other missions, but 
they are presently buried in the classified literature. 

So most of the efforts reported here have concerned 
attitude control as a special discipline, and attitude 
control systems as special major equipment subsystems 
contained in most satellites and space vehicles. Three 
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papers have taken this general point of view in outlining 
the major constituent problem areas of attitude con- 
trol, (16) and (20) as qualitative or summary reviews, 
and (14) as a preliminary quantitative treatment (as 
well as survey) of these areas. The topics discussed in 
those references provide a framework for the literature 
survey below. 


First Principles 


There are three major areas which constitute the 
first principles of the attitude control problem. The 
first is a collection of basic questions and decisions that 
may be lumped together as “foundation problems.”’ 
The second is the quantitative description of attitude 
motions; that is, the formulation of the basic equations. 
The third area is that of performance requirements. 

Perhaps the first basic question to be faced as a 
“foundation problem,” mentioned in (1), (14) and (20), 
is whether attitude control is to be exercised on the 
entire vehicle or individually on a number of its moving 
parts. The answer depends largely on whether the 
vehicle is on a single-purpose or an elaborate multi- 
purpose mission. As (1) and (20) observe, though, some 
control of the whole vehicle is necessary in any case, and 
this is the problem of concern to most of the literature. 
The alternatives of stabilizing several separate parts, 
though mentioned, has not been treated except in a 
work in preparation by the present author. 

Next there is the important question of an attitude 
reference system. Attitude control is not meaningful 
unless it is specified with respect to a definite orientation 
reference. This question is listed in (1), (16) and (20), 
but (14) provides the most complete treatment thus 
far. The latter concerns only a satellite, one of whose 
axes is kept pointed “down” toward the mother planet. 
The relative merits of this orientation and the inertially 
fixed orientation provided by spin stabilization are 
given in (30). A basic reference system for spacecraft 
(as distinct from satellites) is proposed by (31). Al- 
though that paper is more concerned with geometric 
location in the solar system, obviously the frame used 
might be a natural attitude reference. 

To complete the foundation problems, one must 
choose a set of reference axes in the satellite body and 
define suitable attitude deviation angles between these 
and the attitude reference axes. This is done explicitly 
in (14). The center of mass of the body (therefore the 
natural origin of body-fixed axes) wanders because of 
internal mass shifts (3), and the orientation of the 
principal axes of inertia changes within the body for the 
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same reason (10). Both of these effects are developed 
completely, from an analytic point of view, in (17), 
(18) and (22). 

An early step in the quantitative treatment of atti- 
tude control is the development of differential equa- 
tions of motion in which the attitude deviation angles 
are the dependent variables. This is done relative to a 
“natural” satellite attitude reference in (14) for motion 
in three dimensions. For motion in pitch (the plane of 
the orbit) alone, complete equations are developed in 
(28). It has been usual to reduce these to the special 
case of small attitude deviations (14), (21), but the 
cross-coupling of modes and other nonlinear effects for 
larger attitude oscillations has been only a problem 
observed (7), (14), (23), and not yet analysed. 

Several papers (5), (7), and (1) (in the jet control 
case) omit dynamical cross-coupling terms from the 
equations of motion, although these exist even for 
small angles when one uses an attitude reference frame 
that rotates, following a satellite around its orbit. There 
is nothing incorrect in the result if it is realized that 
this omission therefore implies a set of reference axes 
fixed in inertial space, which are natural for some space- 
craft or spin-stabilized satellites, though not for the 
satellite always keeping one reference axis “‘down.” 
This point should be emphasized: although pitch 
motion is independent for small orbits, roll and yaw are 
essentially coupled in the rotating reference frames 
described above. 

An interesting variation of the equations of motion 
is given by Kooy (10). In his treatment, it 1s considered 
that the attitude angles are prescribed functions of 
time and it is desired to find the flywheel speeds (used 
for control; see later) needed to maintain this attitude. 
The flywheel speeds therefore are taken as the depend- 
ent variables. As would be expected, the resulting 
equations look much like those obtained following (14). 
This point of view is particularly appropriate to space- 
craft following a nonperiodic trajectory. It does not 
take account of changes in time of flight (because of 
initial guidance errors, for example) or of perturbing 
torques which generate unexpected attitude motions. 
However, it can be seen with a little thought that this 
approach could be used to determine a nominal base 
state for the control wheels, and that other torque com- 
pensation could be analysed by variation about this 
base state. 

A complete unification and generalization of the 
questions of reference system, attitude deviation 
angles, and quantitative formulation has been at- 
tempted in (23). Equations of motion are developed in 
terms of arbitrarily large attitude deviation angles 
relative to an arbitrary reference frame. For this reason 
it can be said that these general problem areas are now 
in good shape. The most important work still remaining 
is to apply this unified theory to the variety of special 
choices of attitude reference and deviation angles that 
may be of practical interest, and thereby develop the 
basic differential equations for these special cases. 
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Performance requirements of attitude control sys 
tems have not yet received enough attention. The 
nature of the requirements problem is discussed briefly 
in (16) and (20). It is mentioned there that the effect of 
attitude errors often takes the form of a time relation- 
ship. That is, attitude deviations do not represent an 
absolute “loss” in some missions, but merely an ex- 
tended length of time, to gather a given quantity of 
data. But in other cases, such as communication, atti- 
tude deviations may distinctly decrease the overall 
value of the operation. ~~ : 

A more searching examination of missions is highly 
desirable to establish the effects of attitude deviations, 
the criteria which should be used in controlling attitude 
(e.g., control of steady state attitude oscillations vs. 
control of transient attitude peaks), and the explicit 
quantitive requirements on attitude motions. Of course, 
there also are many other requirement factors to be 
satisfied besides accuracy: such as reliability, longevity, 
low power consumption, etc. A good list of general 
qualitative design considerations has been given by 
Haeusserman (7). 


Torques on the Vehicle 


An essential part of the attitude control problem is 
the analysis of the torques acting on the vehicle. Some 
of these will act to perturb the satellite from its desired 
orientation, others to restore it. They arise from a large - 
number of physical effects, many of which are listed 
and discussed in general terms in (14) and (20). 

One class of torques arises from an interaction of the 
satellite with the electromagnetic field through which 
it passes. These are analysed from an elementary point 
of view in (14) for satellites held in a constant attitude 
relative to the earth, but any detailed analysis is very 
sensitive to the exact vehicle configuration. In (32) and 
(33), Vinti develops the torque on spin-stabilized, con- 
ducting spherical satellites in equatorial and non- 
equatorial orbits respectively. He finds not only a 
dissipative component tending to slow the spin, as 
would be expected, but also components normal to the 
spin axis which would tend to cause axis precession. 

The torque source most extensively treated is that 
of the gravitational field. It now is well known that a 
torque exists on an elongated body traveling in a 
satellite orbit. Its quantitative nature has been analysed 
or discussed in (5), (8), (9), (12), (28) with the vehicle 
idealized as one or more pairs of discrete masses ar- 
ranged as “dumbbells” centered on the center of mass. 
This is a permissible representation of the vehicle 
principal moments of inertia as a simple mass distribu- 
tion, but the more general viewpoint of an arbitrarily 
distributed rigid body (of “small” dimensions relative 
to the mother body) is taken in (13). There a generaliza- 
tion is made to an oblate mother body of the MacCul- 
lagh formula for the mutual potential of gravitating 
bodies. On the basis of this result, (15) finds gravitational 
torque components along the principal axes of the body. 
Explicit expressions are obtained for the small-angle 


h se, while the results are left implicit when angles are 
bitrarily large. 

Probably the most important torque of all is the 
i ction torque from internal moving parts, particularly 
pr spacecraft or distant satellite (where the gravita- 
onal torque is small) and for any manned vehicle. The 
Jase of rotating flywheels about the principal axes, or 
, }xes near them, is displayed in (10). A quite general 
eatment of an arbitrary number of internal translating 
nd rotating parts is given in (17), briefly supplemented 
y y (18). A major restriction is that the motion of the 
‘parts is described explicitly by known time functions, a 
\fituation which does not hold when the parts are used 
or control. An extension to the latter case is made in 
22), it being assumed there merely that the translation 
‘hnd rotation of parts is connected through a set of 
‘bontrol equations with the attitude angles (or rates) of 
ithe vehicle structure. In all of these cases, the basic 
‘result is a set of analytical expressions for the torque 
. rom the moving parts. 

| Several other torque sources are discussed as well in 
(14), including the radiation pressure from radiation on 
land from the vehicle, the bombardment of incident 
{particles (ions, molecules, cosmic rays, meteoroids, etc.), 
jand the inertial reaction torque of the coordinate system. 
1The latter is the negative of the torque that must be ap- 
Iplied to the vehicle to make it keep up with its reference 
jcoordinate frame if this frame rotates in inertial 
space with a non-constant angular velocity. (It can be 
shown that the latter is generally the case. See (14).) 

Not many illustrative calculations of torques are 
available, primarily because they are very dependent 
on the exact configuration and orbit. Nevertheless, some 
‘are given in (5), (20), (32) and (33). 

In interpreting the effect of torques, both their 

magnitudes and their time behavior are important. 
Constant torques are particularly troublesome because 
they can saturate the attitude control system (14), (20). 
The gravitational torques depend on the vehicle atti- 
tude rather than time or position, and imply that cer- 
tain vehicle configurations are inherently stable in 
certain orientations. Both (14) and (20) elaborate this 
point and describe how this inherent stability can be 
‘important in the attitude control problem despite the 
low torque levels involved. The free vibrations of such 
an inherently stable vehicle are discussed in (2), but the 
natural frequency of vibrations out of the pitch plane 
is given incorrectly because of the neglect of roll-yaw 
coupling. (It is not difficult to show that no sustained 
free oscillations will continue in any plane except the 
plane of the orbit.) 

Several of the references have considered the torque 
from the gravitational field gradient particularly from 
the point of view of its effect on internal instruments. 
The detection and utilization of the gradient as an atti- 
tude reference direction is emphasized in (3), (8), (9) and 
(12). In fact, most of (12) is devoted to engineering 
discussion of physical devices sensitive to this torque. 

In spite of the fact that much attention has been 
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given to torques on the vehicle, much still remains to 
be done. The gravitational torque is well in hand, as is 
the foundation work for internal moving parts. How- 
ever, explicit torque components need to be worked out 
for a number of specific physical devices involving 
moving parts. Most of the other torques cannot 
be examined in detail in the absence of a specific config- 
uration, but it will be important to lay the foundations 
well enough that the torques can be computed eas- 
ily once a particular vehicle is specified. Another 
worthy effort will be to compile a complete set of design 
principles which will guide the satellite designer in re- 
ducing undesirable and enhancing desirable torques. 


Control Devices 


It is pointed out in (14), (16) and (20) that most 
physical sources of attitude perturbation torques are 
also potential sources of control torques. A general dis- 
cussion of several alternatives is included there. 

Of the non-inertial devices, the possible inherent 
vehicle stability in the gravitational field is the most 
important. However, this importance is due really to 
its operation as a sensor (see later below) rather than 
the actual control torques it provides. 

So it is really the inertial reaction control devices that 
bulk largest in attitude control. Jet control is mentioned 
in a number of references, and is amplified in (1). (How- 
ever, this reference incorrectly omits dynamical cross- 
coupling in the attitude equations of motion which 
occur even when there are no internal rotating parts.) 

On the basis of flexibility and efficiency, reaction fly- 
wheels are very attractive. They are described in (6), 
(10), (14), (21), (22), (26), and (29) among others. 
Reference (22) derives the detailed torque equations 
for identical roll and yaw wheels, (21) for roll and yaw 
wheels that differ. Reference (29) is devoted to pitch 
control by a flywheel. Generally, three wheels are as- 
sumed, one about each principal axis. However, it has 
been observed (1), (10) that because of the roll-yaw 
coupling, two flywheels generally suffice. 

The use of a constant speed pitch flywheel, or a 
pitch wheel which accelerates about a biased value, to 
augment the yaw stability is discussed in (14). Other 
inertial devices have been mentioned (14), (20), (21), 
such as gyrostabilizers and vibration absorbers (equa- 
tions for a type of the latter being given in (22)). Com- 
bined flywheel-jet devices are mentioned in (1) and 
(7). Reference (22) unifies the analytical basis of all 
inertial control devices based on internal moving parts. 

As for status, the analytical foundations for inertial 
methods is well in hand. However, only accelerating 
flywheels, and to a lesser extent jets, have been given 
much physical structure. In general, the time is ripe 
for the design of other physical devices based on these 
principles, with particular attention to other kinds of 
internal moving parts. 

Sensing Devices 


Attitude sensing is discussed from a general point of 
view, i.e., with emphasis on its role, in (14), (16) and 
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(12). In addition, (7) displays a very good tabulation 
of possible attitude sensing methods. One can classify 
sensing devices in various ways: (2) divides them into 
active and passive methods; (20) divides them into 
inertial, line-of-sight, and ambient field. 

The use of gyro-stabilized platforms has been men- 
tioned as one example of inertial methods. This is dis- 
cussed briefly in (29), which points out the principal 
objection to the use of gyros—the drift of the reference. 
It is generally accepted that direct orientation informa- 
tion from free gyros is useful only under two conditions. 
Either the gyros are used strictly as a short-time refer- 
ence, or they are supervised by some independent 
method to correct long-time drift—using the gyro sys- 
tem as a kind of memory, in effect. (See (29)). One 
possible exception occurs when one or more attitude 
modes are already stabilized independently, and then 
the remaining mode can be sensed inertially. A specific 
case given in (20) is the yaw gyro, acting something 
like a gyro compass, when the vertical is stabilized by 
other means. 

And there also is one other major gyro application, 
namely the use of rate gyros to provide attitude damp- 
ing. This is mentioned in a number of places, including 
(1) and (20). One problem associated with rate gyros 
is their inability to sense merely the deviation rates 
from the attitude reference system. They also pick up 
the rotation of the reference frame itself, and this must 
somehow be biased out before the signal indicates devi- 
ations only. The latter question deserves further atten- 
tion from the system design point of view. 

The line-of-sight method for attitude sensing most 
often is horizon scanning, appropriate for satellites or 
any spacecraft relatively near a planet. Principles and 
some discussion of errors of one device of this kind is 
the subject of (19). Both (28) and (29) state the impor- 
tance of the method, but analytical details (beyond 
(19) are sparse, and detailed descriptions of physical 
devices do not seem to have been published. 

Other line of sight possibilities have been mentioned. 
For example, (11) gives radio techniques as a possible 
method of attitude indication, although the use was 
not explored. The present author believes that this 
problem deserves considerably more attention. Among 
the miscellaneous methods are direct passive sun 
sights (to give two attitude angles) and photographs of 
sun and planets against fixed stars (31). Although the 
latter reference does refer briefly to attitude control, it 
is concerned almost exclusively with the determination 
of spacecraft location in the solar system. It would not 
be included among the present attitude control refer- 
ences at all except that the title can be easily misinter- 
preted. By including it, together with these comments, 
the reader is explicitly warned that it is not much con- 
cerned with attitude control. 

The methods based on the sensing of ambient fields 
have been devoted almost exclusively to the sensing of 
the gravitational field gradient (mentioned previously 
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in connection with torques on the vehicle). An excep 
tion is the method of Stuhlinger (29) utilizing cos mit 
rays incident on the vehicle. Although this is an i 
genious method, there is a question whether, on th 
basis of existing data, one can demonstrate adequat; 
ray collimation for seniude sensing purposes. A seriou 
evaluation of the method probably must await furthes 
experimental evidence. | 

The true role of the gravitational field gradient acting 
on the entire vehicle, namely attitude-indication 09 
sensing rather than the realization of significant contro) 
torques, is emphasized in (14). Most other discussions 
(3), (8), (9), (12) are more concerned with the action o HI 
internal devices which then might be used to sense the 
direction of the field gradient. In (9), a sensing device is 
proposed consisting of two internal pendulums with 
damping between them. A single pendulum is described 
in considerable engineering detail by Oberth in (12). 
The treatment of the pendulum occupies the bulk of the 
paper; in fact, it constitutes the only real engineering 
of a physical device discussed in any of the references. 
Some difference of opinion exists as to the degree to: 
which this particular method will be useful, but the 
present author can only hope that the same careful) 
engineering attention will be given to other possibilities: 
as Dr. Oberth has given to this one. | 


Attitude Control Systems 


The real piéce de résistance of attitude control, of 
course, is the complete control system. It constitutes a 
control or actuator device, generally an attitude or 
attitude rate sensing device, and a control filter or 
computer. The fundamental central system problem is 
the synthesis of this control filter, as described in (21). 
Although a number of system types have been men- 
tioned in the literature, not many details of the syn- 
thesis-or performance analysis have been presented. 

It is appropriate just to list some of the types pro- 
posed: 

1. spin stabilization 

This is the method used for the Vanguard and 

Explorer satellites. Some of its advantages and 

disadvantages are listed in (29), but the most 

complete analysis of a system (a lunar instrument 
carrier) using spin stabilization has been developed 

by Buchheim (4). 

. reaction flywheels 

These have been mentioned as sole devices by 

Esnault-Pelterie (6) and Oberth (12). 

3. reaction flywheels plus an augmentation system 
The augmentation system can prevent torque or 
angular momentum saturation, usually taking the 
form of a set of reaction jets. This method is 
favored by Angle (1) and Haeusserman (7). The 
latter is one of the few current works containing a 
number of details on the system synthesis and 
performance. It also is generally assumed that 
damping is present in the system. 


bo 


4. reaction flywheels with damping 
Oberth (12) has suggested both magnetic and 

_ fluid internal damping to accompany a flywheel 

. system. 

5. inherent vehicle stability 

_ This possibility is considered in some detail in (28), 
which is another of the few references analysing 
performance in any detail. The case considered is 
one-dimensional (pitch) control, and the treat- 
ment utilizes computer solutions. 

6. inherent stability plus “natural” internal damping 
Reference (2) proposes this as a very simple 
method, but points out the disadvantage that 
attitude errors damp out very slowly. 

7. inherent stability plus reaction flywheels plus rate 

gyros. 

This system is probably the most sophisticated 
and flexible currently proposed. It is described in 
(1), and general system equations are given in 
(22). In addition, (21) presents some actual results 
obtained by computer solution, although it 
passes over the procedure used to choose system 
parameters. 

In terms of general control system characteristics 
ather than specific methods, the following points have 
een made by various authors. First, if the vehicle con- 
ains a net internal angular momentum of moving 
arts, it is important that this have no projection on 
he roll or yaw axes (14). At the same time, this net 
ternal angular momentum (perhaps deliberately 
creased) can be used to enhance the vehicle’s yaw 
tability (1), (14). Second, the possible importance of 
onlinear control with a dead zone about the reference 
ttitude, thereby saving control power, is stressed in 
7). Third, the inherent gravitational stability is dis- 
ussed in terms of vehicle configuration in (14), and 
ome implications are noted. 

As mentioned above, very little attention has been 
iven to the details of rational control system synthesis. 
methodology is outlined and discussed in (14), (16) 
nd (20); moreover, (21) goes just a little deeper into 
he synthesis problem for a specific system type (num- 
yer 7 above) and presents some results, but it skips 
nany of the intermediate steps. 

There are three areas where work is especially desir- 
ible, and these are among the most important in the 
vhole subject of attitude control. One is to devise other 
nechanization schemes, other combinations of devices, 
or control. The second is to consider more seriously 
ind much more thoroughly the synthesis problem as a 
sroblem of servomechanism design. And the third is to 
nalyse the resulting system performance with special 
ttention to the effects of instrument errors. 


Jther ‘Topics 


There are other important questions which are not 
neompassed by the preceding outline. One of these is 
he role of man in attitude control. A provocative dis- 


cussion of some of his pontentialities in this area has 
been given by Senders (24). It is felt that this is an 
area ripe for further investigation. 

Another is the general area of design principles to 
ease the attitude control problem as much as possible. 
Some of these principles could be listed out of hand, 
others could be developed with further investigation. 
But in this paper one major point in this connection is 
emphasized, quoted from (20). ““From the point of view 
of attitude control, a satellite or space vehicle is not a 
structure, but a precision instrument!” 

It is hoped that this summary of the attitude control 
literature and the major problems remaining will 
stimulate more activity on this important and inter- 
esting subject. 


References 


1. Anate, E. E., “Attitude Control Techniques,” (portion 
of “Attitude Control: A Panel Discussion’’), Naviga- 
tion, 6 (1958), 66-71. 

2. Baker, R. M. L., Jr., “Passive Stability of a Satellite 
Vehicle,” (portion of “Attitude Control: A Panel Dis- 
cussion”), Navigation, 6 (1958), 64-65. 

3. Bont, A., “Artificial Satellite, Unification and Mechanics 
(Sidar-Mechanies),”’ Astronautica Acta, 1 *3 (1955), 
120-136. 

4. BucuHerm, R. W., “Lunar Instrument Carrier-Attitude 
Stabilization,’ RAND report RM-1730, 4 June 1956. 

5. Davis, Wo. R., “Determination of a Unique Attitude for 
an Earth Satellite,” presented at 4th Ann. Mtg. Am. 
Astron. Soc., N. Y., 30 Jan. 758, AAS preprint 57-10. 

6. Esnautt-PeLtTerin, Roper, ‘“‘L’Astronautique,” Im- 
primerie A, Lahure, Paris, 1930. 

7. HanussprRMANN, W., “An Attitude Control System for 
Space Vehicles,’ Presented at the ARS Semi-Annual 
Meeting, Los Angeles, June 9-12, 1958; ARS Reprint 
642-58. 

8. Kuemprrer, W. B., AND Baker, R. M. L., Jr., “Satellite 
Librations,” Astronautica Acta, 3 (1957), 16-27. 

9. Kuemprrer, W. B., “Satellite Librations,” (portion of 
“Attitude Control: A Panel Discussion’), Navigation, 
6 (1958), 62-63. 

10. Kooy, J. M. J., “On Plotting Small Thrust Space Ship 
Orbits,” in Probleme der Weltraumforschung, Proc. 4th 
Int. Congr. Astronautics, Zurich, 1953, 107-114. 

11. Merren, R., ‘Uber einige hochfrequenztechnische 
Probleme der Weltraumfahrt,” in Probleme aus der 
Astronautischen Grundlagenforschung, edited by H. H. 
K6élle (Vortrige III Internationalen Astronautischen 
Kongress, Stuttgart, 1-6 Sept 1952), 289-244. 

12. Osertu, H., “A Precise Attitude Control for Artificial 
Satellites,” in Vistas in Astronautics, (International 
Series of Monographs on Aeronautical Sciences and 
Space Flight, Division VII, Astronautics Div., Vol. 1), 
Pergamon Press, London, pp. 217-255. 

13. Roperson, R. E. anp TatisrcHerr, D., “The Potential 
Energy of a Small Rigid Body in the Gravitational 
Field of an Oblate Spheroid,” J. Frank. Inst., 262 
(1956), 209-214. 

14. Roperson, R. E., “Attitude Control of a Satellite 
Vehicle—An Outline of the Problems”, Proc. VIII 
Intern. Astronaut. Federation Congr. (Barcelona, 1957 
Springer-Verlag, Wien, 1958, 317-339. 

15. Roperson, R. E., “Gravitational Torque on a Satellite 
Vehicle,” J. Frank. Inst., 265 (1958), 13-22. 

16. Roperson, R. E., “Problems and Principles,” (portion 


The Journal of the Astronautical Sciences 29 


of “Attitude Control: A Panel Discussion”), Naviga- 
tion, 6 (1958), 59-61. 

17. Rosurson, R. E., “Torques on a Satellite Vehicle from 
Internal Moving Parts”, J. Appl. Mech., 25 (1958), 
196-200. 

18. Rosurson, R. E., “Torques on a Satellite Vehicle from 
Internal Moving Parts (Supplement)”, J. Appl. Mech., 
25 (1958), 287-288. 

19. Rosurson, R. E., “Optical Determination of Orientation 
and Position Near a Planet,” Jet Propulsion, 28 (1958), 
747-750. 

20. Roperson, R. E., “Where Do We Stand on Attitude 
Control?”, Aviation Age, Research and Development 
echnical Handbook, June 1958, B-5 to B-10. 

21. Rogerson, R. E., ‘Inertial Control of Satellite Attitude” 
Elect. Eng., TT (1958), 1086-1088. 

22. Roperson, R. E., “Principles of Inertial Control of 
Satellite Attitude”, presented at the IX Intern. 
Astronaut. Federation Congr., Amsterdam, Sept. 
1958. 

23. Roperson, R. E., “A Unified Analytical Description of 
Satellite Attitude Motions”, presented at the Annual 
Meeting Am. Rocket Soc., New York, Nov. 1958. 

24. Senpers, J. W., “Human Factors in the Design of 
Guidance and Control Systems for Missiles and Space 
Vehicles,’ ARDC Symposium on Guidance of Ballistic 
Missiles and Space Vehicles, WADC, 11-13 March 
1958. 

25. Sincmr, 8S. F., “Application and Design Characteristics 
of Minimum Satellites,’ ARS reprint 278-55 (Nov. 
1955). 

26. Smivu, F. A., “The Satellite Telescope,” J. British Inter- 
plan. Soc., 16 (1958), 361-367. 

27. Srruuine, K. R., “Earth Scanning Techniques for a 
Small Orbital Rocket Vehicle,” Probleme der Weltraum- 
forschung, 4th Astronautical Congress, Zurich, 1958, 
63-70. 

28. Stocker, T. A. J. AND VacHINo, R. F., “The Two- 
Dimensional Librations of a Dumbbell-Shaped Satel- 
lite in a Uniform Gravitational Field,” MS thesis 
School of Engineering, Air Force Inst. of Technology 
Air University, WPAFB, March, 1958. 

29. STUHLINGER, E., “Instrumentation Problems of Un- 
manned Satellites,” ARS Repr. 806-56 (June, 1956). 

30. SrunLineer, E., “Control and Power Supply Problems 
of Instrumented Satellites,” Jet Propulsion, 26, *5 
(May, 1956), p. 364. 

31. Vertreet, M., “Orientation in Space,” J. Brit. Inter- 
plan. Soc., 15 (1956), 324-338. 

32. Vintt, J. P., “Theory of the Spin of a Conducting Satel- 
lite in the Magnetic Field of the Earth,’’ BRL Report 
1020, Aberdeen Proving Ground, Maryland, July, 
1957. 

. Vintt, J. P., “Theory of the Spin of a Conducting Satel- 
lite in Non-Equatorial Orbits,’ BLR Report 1031, 
Aberdeen Proving Ground, Maryland, October, 1957. 


Addendum 


Farrior (84) recently has discussed in more detail 
some of the reasons why attitude control is needed for 
many missions. The discussion is qualitative without 
angular accuracy requirements being given. However, 
it was the consensus of a panel on “Space Vehicles, 
Satellites and Missiles” (for which (21), (34) and (35) 
are the prepared comments) that accuracy requirements 
cannot be given in general terms. They are strongly 
mission-oriented and may vary from a small fraction of 
a second of are to 20° or 30°—perhaps more. 


we 
we 
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An additional perturbation (and potential cont 
torque source is treated in (35), namely air drag. i 
| 

C 


the gravitational gradient effect discussed earlier ab 
air drag causes a torque that depends on the orientati 
and geometry of the satellite. Particular attention 
given in (35) to identifying the altitude regime whe 
air drag torques dominate gravitational gradier 
torques and conversely. Because the positions of stab 
equilibrium differ for these torques, their combi c 
effect in the transition region where they have con 
parable magnitudes is to shift the equilibrium positic 
to an interesting and important phenomenon f 
relatively low satellites. 

One additional discussion of the attitude contr 
problem for a specific mission, a lunar circumnavig 
tion, has appeared (36). It is proposed there that th 
attitude be stabilized during surveillance of the moo: 
using an astrotracker as a primary attitude informatica 
input and reaction flywheels as actuators. Special a7 
tention is given to system reliability, a refreshin 
novelty. The most critical part of the system from th. 
point of view was the postulated digital computer, bu 
other computer reliability estimates make the pictur 
much less bleak. The reliabilities of several alternatiy 
systems are estimated and compared, some of whie 
give acceptable (if not attractive) probabilities c 
success. | 
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Logical Approaches to Problems of 
Space Flight’ 


Commander George W. Hoover, USN (Ret.) 


Office of Naval Research, Washington, D. C.+ 


Introduction 


The American Astronautical Society has by its own 
charter placed itself in a unique and challenging posi- 
tion. By implication, the members have decided to con- 
tribute their joint capabilities toward solving the 
problems of Space Travel through the development of 
the astronautical sciences. A commendable objective— 
but an extremely difficult task.*** 

First of all when we refer to space travel we are talk- 
ing about a man-machine system. Such a system di- 
vides itself into a “‘brain’’ and a ‘‘muscle’’. 

The ‘‘muscle”’ can be broken down into: (1) structure 
and materials, (2) aerodynamics and ballistics, (3) 
propulsion. 

The “brain” can be divided into sensors, computa- 
tion, display and mechanism amplifiers, displays, and 
‘controls. 

In addition to the above we have the problems of: 
(1) communication, (2) aero-medicine, (8) geophysical 
and astrophysical. 


Major Requirements 


Now that the problem is divided into specific areas we 
can establish the requirements for each. This is not an 
easy task but is absolutely essential if we are to deter- 
mine the proper direction in which to look for the com- 
pletely adequate answer. In the brain portion one pri- 
mary question to be answered is simply, ‘What is it 
the pilot needs to have displayed?” 

There are actually two major requirements. These 
are, position in space; and position astrophysical. The 
first is observed normally with respect to a vertical 
plane and the second with respect to a horizontal 
plane. The first is primarily orientation, the second 
situation display. 

Each of the two areas must be further divided into 
three types of display: (1) Orientation, which tells the 


* Editor’s note: Although the Journal of Astronautical 
Sciences primarily emphasizes quantitative treatments, it 
also is open to sound discussions of the foundations or phil- 
osophical aspects of particular problem areas. These remarks 
of Cdr. Hoover are of this nature. They were presented at 
the Summer Space Flight Meeting of the American Astro- 
nautical Society, September 13, 1957, New York, N.Y. before 
the establishment of the first satellite, but they are no less 
timely today. Although some portions of the original address 
unfortunately have had to be deleted because of space limi- 
tations, the essential points of Cdr. Hoover’s address have not 
been changed. 

+ Now Director of ‘Technical Planning, Benson-Lehner 
Corp., Los Angeles, Calif. 


man what he is doing; (2) Director, which tells the man 
what he should be doing; (3) Quantitative, which tells 
the man how he is doing. 

In order to determine the details of these three 
types of display it is necessary to make a complete 
analysis of the Information Requirements. To deter- 
mine these requirements the operators must be inter- 
rogated to establish the fundamental information 
necessary to accomplish the task at hand. Not an in- 
terrogation of the individual’s opinion, but rather an 
assessment of the essential information required with 
which the task can be accomplished without interpola- 
tion or mental computation. For example, during an 
orbit interception the pilot must have range informa- 
tion. Engineers realizing this would provide a means for 
determining range and then proceed to present it in 
miles and yards. However, when prospective space 
pilots are asked the question, ‘‘Why do you want 
range information”’, they answer by stating they need 
range data in order to know when to reduce thrust, 
when to reverse thrust, and when to stop thrust. When 
denotes time, not distance. In other words, fundamen- 
tally, range must be indicated as a director type of 
presentation in order to eliminate the necessity of the 
pilot remembering at which range to reduce, reverse, 
and stop thrust. When the questioning reaches a point 
where the operator states that without X data he cannot 
carry out the task, then X data is the fundamental in- 
formation requirement. *** 

By using the visual world as our yardstick it is rela- 
tively simple to determine what it is in the visual world 
that causes us to react as we do throughout various. 
phases of flight. If we resort to invention we can only 
compare one idea with another, both of which may 
be fundamentally wrong. On the other hand, if we use 
the visual world as one of the elements of our compari- 
son, which we know incidentally is adequate, we at 
least can establish an equally adequate display. This 
comparison then is made against a natural, rather than 
a man-made model. If we can create a synthetic display 
comparable, (not a duplicate) to the visual world, then 
we know axiomatically that the display will be ade- 
quate. Having once established this we can then pro- 
ceed to bring about innovations to take care of the de- 
ficiencies in the visual world which are apparent during 
the special modes of space flight. 

Such an approach to determining a display is merely 
resorting to the concept of “doing what comes natu- 
rally’. We have no difficulty in walking in clear 
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weather nor in driving a car in clear weather. By the 
same token we fly well in clear weather. This then 
should be the display for which we must create a corre- 
late. The proof here for a display is by axiom rather 
than evolution because if it is true that we orient our- 
selves by our perception of the visual world, it follows 
that our display must be adequate if it reproduces the 
same cues which are apparent in the visual world.*** 

If conventional methods of conducting research and 
development in this field are carried out the answers 
may be forthcoming, but not for a long time, and not 
without a tremendous expenditure of money. There 
are many reasons why conventional approaches to the 
problems of research and development are inadequate 
but perhaps the most predominant is the influence of 
inhibited thinking.*** 

Perhaps the most important cause of inhibited think- 
ing, however, is due to the methodology employed in 
carrying out most research and development. Too much 
reliance is placed on intuition and invention. Too often 
the expert opinion is really not much more than per- 
sonal opinion, or where a group is involved, a concensus 
of opinion. Because an engineer or scientist has worked 
in a particular field for a long time does not in itself 
qualify him to stand unquestioned. Too often decisions 
are made as a result of past experience which may have 
been excellent for previous circumstances but will not 
hold up under new requirements. Such opinions are not 
always backed up by facts. The expression, ‘‘Don’t 
confuse me with the facts, I’ve already made up my 
mind’’, is not always made in jest. 

Programs and projects are started too often without 
making a complete statement of the problem which 
results in only a partial solution which is justified by 
calling it a compromise. 

An approach which permits inhibited thinking can 
only lead to short term answers with little or no future, 
ideas resulting in inventions requiring lengthy develop- 
ment and test programs, opinions often leading to blind 
alleys, and small improvements at a very large cost. 


Approaches to Problem 


Actually we have had the right approach with us 
since 1620 when Sir Francis Bacon in his Novum 
Organum made the following statement: 

“It is idle to expect any great advancement in 
science from the superinducing and ingrafting of new 
things upon old. We must begin anew from the very 
foundation, unless we would revolve forever in a circle 
with mean and contemptible progress.” 

In other words, we must stop modifying modifica- 
tions. We must stop accepting compromises or partial 
solutions. We must think uninhibitedly and seek an 
adequate solution by stating the problem in its funda- 
mental terms. In order to think in this manner we must 
treat each problem as a completely new one and start 
by wiping the slate clear. 

Such an approach calls for the user to seek not a 
partial solution, but rather an ultimate solution. Many 
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engineers today consider such tactics as dreaming. 
Engineers seem to prefer to be what they call objective 
in their thinking, practical, with their heads not in the 
clouds, but with their feet solidly on the ground. This 
is certainly a very healthy attitude, providing their feet 
don’t get stuck too deeply in the mud.*** 

From the information requirements and the display 
requirements we can determine that the following 
sensors are necessary: (1) inertia, (2) force, (3) heat, 
(4) quantity, (5) geometry, (6) electromagnetic radi- 
ation. ; : 

All of the equations related to the flight of any space — 
craft can be solved by using these six sensors with — 
respect to time. Navigation is a function of velocity 
which is a function of acceleration. Fuel management. 
is a function of fuel available with respect to time. 
Position, obstacles, and weather are some of the infor-— 
mation requiring the use of some form of electro- 
magnetic radiation. Here again these sensors are. . 
fundamental because they represent the basic param- | 
eters relevant to the motion of any craft in space. They 
are all variables in the equations of navigation and | 
orientation and therefore all modes of flight. 

Since information requirements cannot be deter-— 
mined, these are replaced by operational or functional } 
requirements. In other words, what is the machine re- | 
quired to do? Following this we must determine the | 
environmental requirements or—in what environment 
must the machine operate. Next we determine the struc- 
tural or configuration requirements. In this instance 
we establish what general shape or design the machine 
must conform to in order to meet the operational and 
environmental requirements. And finally we must de- 
termine the material requirements in order to meet the 
previous requirements. In other words, what must the 
materials be? 


Concluding Remarks 


What I am suggesting is not an easy thing to do be- 
cause it necessitates divorcing from your minds every- 
thing you have learned except the fundamental laws of 
physics during the determination of the requirements. 
The worse pitfall during this process is—falling in love 
with new possibilities which will arise—and thus dis- 
continuing the seeking out of fundamental. 

What I have suggested here may be quite confusing 
but I assure you that, based on experience, such an ap- 
proach will work. When the fundamental requirements 
of any problem are established, the solution to the prob- 
lem becomes apparent. If you will take the time to state 
the problem at hand in its fundamental terms, if you will 
stop inventing and seek a completely adequate solution, 
if you will look at the total problem and work asa team: 
then you will be thinking uninhibitedly and progress will 
be made not in small increments, but in a continuous 
series of major breakthroughs. 

In simple words: state the problem completely before 
you try to solve it and space travel will become a reality 
rather than a discussion. 


Format of Technical Papers 


for AAS Meetings 


At the suggestion of various members of the 
Society the Technical Papers Committee Chair- 
man, Dr. Horace Jacobs, has prepared a sample 
sheet of the basic requirements for technical 
papers to be submitted for possible inclusion in 
various AAS meetings. This is reproduced at this 
time for wide distribution. Abstracts of such com- 
munications should be typed single-spaced in a 
text width of 4.5 inches centered below the title 
and author(s) of the paper. The abstract can be 
brief and should not exceed 300 words in length. 
1. Introduction. A technical paper for publication 
in the AAS Proceedings should be typed on a 
good quality bond, 8} x 11 inches. Except for 
the abstract, page width is 63 inches and page 
length is 9 inches. Copy should be typed in 
standard elite, preferably on an electric type- 
writer, and should be suitable for offset. Cor- 
rections should be stripped in rather than 
erased. Text is single-space with double-space 
between headings, equations, items, and para- 
graphs. Major heads are all caps and flush left. 

2. Other Main Headings. The paper can be divided 
into principal sections as appropriate. Headings 
or paragraphs are not numbered. Secondary 
headings are in caps and lower case; they are 
flush left and underlined. 

a. Equations. Equations are as follows: 


(a?mo/ ats?) (a?mo/as2?) = 0 (1) 


Slashes should be used to separate numera- 
tor and denominator if it is conveniently 
possible. Equations are referenced in the 
text as equation (1), (2), ete. 

b. Symbols. Slashes are also preferred to 
separate numerator and denominator in the 
text: m:; = mo exp (—v;/c). Where super- 
scripts and subscripts make single-space 
typing impractical, the lines where they 
occur can be typed in space-and-a-half 
or in double-space. 

c. Greek Symbols. If tertiary headings are 
necessary they are flush left, upper and lower 
case, and underlined; text is run-in. Greek 
symbols can be prepared in one of the follow- 
ing ways: 

(1) Typed by varitype 


(2) Typed by a mathematical typewriter 

(3) Hand-drawn. If this method is used the 
following should be noted: Use black or 
red ink, Draw symbols neatly. 
Above is an example of format for enumer- 
ations. They are indented 5 spaces. 

d. Tables. Table headings are centered above 

the table and have the following format: 


Table I 
THE HEADING IN ALL CAPS 


e. Figures. Line drawings suitable for reproduc- 
tion may be inserted by the author at 
appropriate places in the text. If this is done, 
figures should be reduced consistent with 
clarity and space conservation. Captions 
are centered below the figure and have the 
following format: 


Fig. 1. The Caption of the Figure 


References to figures are made in the text 
by the designation Fig. 1, Fig. 2, etc. 

If line drawings are not inserted in the text 
they should be provided as original artwork 
or photostats suitable for reproduction. 
Glossy black and white prints should be 
provided for half-tone photographs. Each 
figure should be numbered and a caption 
provided. 

f. References and Footnotes. References listed 
at the end of the paper are indicated in the 
text by a superscript Arabic number’. 
Footnotes are called out in the text by super- 
script symbols *, **, ft, {, etc. Standard 
practices footnotes and references are 
acceptable. 


3. Conclusion. The AAS would appreciate it if 


each author would provide the Society with 
at least 100 preprints of his technical paper 
before the beginning of the session. These 
preprints should be presented in the format 
indicated above. However, illustrations 
might be interspersed in the text or placed 
at the end of the paper. 
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